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Abstract 

We investigate the nonlinear dynamics of two coupled annular Bose-Einstein condensates (BECs). 
For certain values of the coupling strength the nonrotating ground state is unstable with respect to 
fluctuations in the higher angular momentum modes. The two branched Bogoliubov spectrum ex- 
hibits distinct regions of instability enabling one to selectively occupy certain angular momentum 
modes. For sufficiently long evolution times angular momentum Josephson oscillations sponta- 
neously appear, breaking the initial chiral symmetry of the BECs. 

PACS numbers: 03.75.Lm, 03.75.Mn, 03.75.Kk 



1 



One of the most famous paradigms of quantum physics is the existence of Josephson 
oscillations. They were first predicted for superconductors separated by an insulating layer 
Later they have been observed in superfluid '^He ^ and gaseous BECs 0, 4|. Next 
to the oscillations of charge and/or particles between two modes these systems can exhibit 
highly nonlinear dynamics with sometimes surprising behavior. In this letter we study one- 
dimensional (ID) BECs confined in two ring-shaped traps which are sufficiently close to 
each other to allow tunneling through the barrier between them. We demonstrate that the 
stationary state, in which only the zero angular mode is occupied by the BECs, becomes 
unstable for certain values of the coupling strength. This allows for a selective occupation of 
specific angular momentum modes. For short propagation times the angular momentum in 
each ring is conserved. For longer interaction times, however, angular momentum Josephson 
oscillations appear. This novel type of Josephson oscillations spontaneously breaks the initial 
chiral symmetry of the individual BECs. 

In the mean-field description the evolution of a dilute gas of identical interacting Bosons 
under the infiuence of the trapping potential V{r) is governed by the Gross-Pitaevskii equa- 
tion (GPE), which in cylindrical coordinates reads 

= - 9! + ^) + ^(r)]^ + 9\^\'^ (1) 

where M is the atomic mass, g the nonlinear coupling constant, \1/ = \l/(r) the bosonic mean 
field and = —h^d"^ the z-component of the angular momentum operator. The system 
consist of two BECs in parallel ring-shaped traps encircling the 2;-axis. The positions of the 
upper and lower ring are ±2:0, respectively. Correspondingly, the trapping potential takes 
the form V{t) = Vp{p — po) + Vz{z, zq). The first term provides harmonic radial confinement 
centered at p = po, and Vz{z, zq) creates a symmetric double well potential with its minima 
a.t z = ±zo. Both BECs reside in the radial ground state \I/p(p) of Vp{p — po). Vertically 
they occupy the harmonic ground states $(2; ± ^o) which are localized in the upper and 
the lower well of Vz{z, zq), respectively. The total wave function of the system can then be 
written as \E'(r) = \I'p(p) [^{z — Zq)xu{<P) + ^{z + ZQ)xd{<P)] where the indices u and d refer 
to the upper and lower ring. After inserting \E'(r) into the GPE we obtain the two coupled 
equations 

idrXu/d = -dlxu/d + i^Xdiu + l\Xu/d?Xu/d (2) 
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Here we have introduced the scaled time r = JMW^^ coupling k = 

—R^j dz^{z + 2:0) [d], — ^V{z)'j $(2; — 20) and the interatomic interaction parameter 
7 = ^^^2^1^ / dz^'^iz) with R-^ = J rfpp-2|^p(p)|2 [l5|. We assume that the external 
confinement allows an effective ID treatment of the BECs. The atom-atom interaction 
can then be described by the ID coupling constant gn) = |6(- Here a is the three- 

dimensional s-wave scattering length and ap the harmonic oscillator length of the radial 
ground state. Equations, similar to eq.([2]) arise in the context of two coupled elongated 
condensates Q] • Ring traps additionally allow the existence of stationary currents. 

In the angular momentum mode representation the azimuthal wave function of the indi- 
vidual BECs can be written according to Xu/d = (271)^^^^ exp(i du/d) Xlm «m'^exp(im0) with 
6u/d being the phase of the wave function in the respective annulus. The coefficients am are 
normalized such that jam'^p = Nm'^ corresponds to the number of particles residing in the 
m-th angular momentum mode. Hence J dcp = N'^^'^ corresponds to the total number 

of particles in each of the two annuli. Inserting the above expression for Xu/d into eq.([2]) we 
find the system of coupled equations 

zdra:(' = rr?at' + ^^id^t + f J] oi'il'ot^'' o^t' (3) 

/vr ^ — ^ 

nn' 

with K^jfi = K*^i^ = K e*'-^'*"^"-'. We now seek a stationary solution of this system for which 
in both of the annuli solely the m = mode is occupied. This exists only for equal number 
of particles, i.e. A^^q = Nq = Nq, and equal coupling = Hd = i^- We then find the two 
solutions 

a- = = v^e^(^±'^)^+^^ , afi^, = (4) 

with some arbitrary phase 6 and e = being the nonlinear energy due to the interatomic 
interaction. Hence the total two-dimensional wave function becomes either a symmetric or 
an anti-symmetric superposition of the axial ground states of the two annuli: 



^±(r) = 0Voe^^±^+^^^,(p) [$(2 + zo) ± $(2 - zo)] 

with fi± = 6 ± n being the chemical potential. 

In order to investigate the stability of these states with respect to fiuctuations in modes 



with m 7^ we make the ansatz a^^Q = e 



Inserting this together 
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with eq. (jlj) into the eq. ([3]) yields after hnearization in the Um and Vm the eigenvalue equation 

cou^;^' = (m^ + /ii)^:;/'^ + ee^2%!{^ + nufi^ (5) 

We then find the excitation spectrum consisting of the two branches 

u)± = a/ (m^ + e — K± — e'^. (6) 

This result is obtained for either chemical potential fi±. The branch corresponds to 
the well-known Bogoliubov spectrum Q] of a uniform BEC but with integer m. In this 
letter we consider only repulsive interatomic interaction, i.e. e > 0. Thus ct;+ is always a 
real number and the coupled condensates are stable against fiuctuations in this excitation 
branch. Note that is independent of the coupling strength k. The cu. branch on the 
other hand depends on k, and result from the interaction among the coupled condensates. Its 
frequencies are either real or imaginary but never complex. For zero coupling U- is identical 
with u^. Figure [1] depicts for two values of the nonlinear energy {e = 1.0 and e = 2.0) 
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FIG. 1: (Color online) Branch of the Bogoliubov spectrum of the coupled condensates plotted 
against the coupling strength k for two values of the nonlinear energy e (blue/solid: real part, 
red/dashed: imaginary part). In panel a the instable regions are labelled with the respective m 
number. Due to the symmetry of eq.([5]) the spectrum is symmetric with respect to a sign change 
of m. 

the spectrum uj_ as a function of the coupling k. The blue/solid curves represent the real 
part and the red/dashed curves the imaginary part of u^. First, let us consider relatively 
low energies {e < 2) (fig.[T^). Starting from k = 0, we notice that is positive and real. 
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and that it decreases monotonously to = at k = 1/2. After this, one enters a region 
where uj- is imaginary with Im(co'_) > 0. The system is then unstable under fluctuations in 
the m = ±1 modes which grow at a rate of F = 2Im(co'_) [?, 9]. As k increases regions of 
stability and instability follow one another. The latter are defined by m? /2 < n < + 
For any given e the maximum growth rate F^ax = 25 is a universal quantity for all modes 
which is independent of m and is established at the coupling strengths k = \[Tn? + e\. 
Note that for sufficiently large e the unstable regions of two adjacent m-modes may even 
overlap thus eliminating the stable region in between. This can be seen in fig.[T)D, where 
for a nonlinear energy e = 2.0 part of the regions of instability for the m = 1 and m = 2 
modes overlap. Since the coupling k is a function of the trapping potential, and as such 
experimentally accessible, the instability of the modes can be used to selectively affect one or 
more m-modes of the rings. For example, at e = 1.0 each mode is "individually addressable" 
through an appropriate choice of k (see fig.[T^). At e = 2.0 and k, = 2.1, on the other hand, 
both the 771 = 1 and m = 2 are unstable with respect to fluctuations. 
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FIG. 2: (Color online) Evolution of the occupation number normalized to the total number 
of particles for the nonlinear energy e = 2.0. For k = 1.6 (panel a) we observe an exponential 
population increase in the m = ±1 modes. Population of the m = ±2 mode is also visible at 
later times but is suppressed by more than two orders of magnitude. For k = 3.2 (panel b), the 
accumulation in population of the m = ±1 modes is suppressed and only the exponential growth 
of the TTi = ±2 modes is visible. 



The growth of the instable modes eventually leads to a break down of the linearized eqs.(l5]) 
due to the interaction between higher lying m-modes. We therefore return to eq.(l3]) and 
integrate them numerically. We choose Oq^'^ = ^JNq + 6u/d as initial condition, i.e. almost 
equal number of atoms in the m = mode of both rings, whilst allowing for experimentally 
unavoidable particle number fluctuations of the order of 6u,d = Adding these 



fluctuations has only minor influence on the numerical results. Following Saito et al. [9|, we 
introduce a small seed in the lowest few angular momentum modes (up to m = ±5) with 
a magnitude of 10~^ x \/Nq. Again, such a fluctuation is experimentally inevitable. We 
truncate the set of coupled eqs.(l3]) at the angular momentum mode m = ±15, well above the 
highest contributing mode. We verified the quality of the propagation by monitoring energy, 
norm, and angular momentum conservation. Since both annuli have a slight population 
difference the nonlinear energy is calculated according to e = with A^tot = Nu + Nd- 
An example of the numerical propagation can be seen in fig. [21 We show the occupation 
A!^ for m = 0,±1,±2 at the nonlinear energy e = 2.0 (see also fig. [lb) for two different 
coupling strengths k. For n = 1.6 only the m = ±1 modes are unstable. We observe for 
early times (r < 10) the predicted exponential increase in population with a rate of F ^ 4. 
At later times r > 9 the population also in the m = ±2 modes increases slightly. This 
cannot be described by the linearized eqs.([5]). However, in the time window considered here, 
the population of the m = ±1 modes is more than two orders of magnitude larger than one 
of the m = ±2 modes. In the lower panel we present the same plot for k, = 3.2. Here we 
observe no population growth within the m = ±1 modes but N!^2 grows at a rate of F = 3.9. 
This clearly demonstrates the possibility of a selective angular momentum mode excitation 
by tuning k. 

Let us now turn to the angular momentum of the two BECs. The Lz expectation value of 
uj_ branch is (i^")™ = (L^)"^ = which implies an equal population of states with opposite 
m. Conversely, for modes of the branch uj+ we find {LD^ = {Lf)_^ = y ^"^2"^^ —- However, 
according to eq.(l6]) the value of uj+ is always real, which results in a growth rate of F = 0. 
In fig. [3] we present the angular momentum per particle and the relative particle number 
difference between the coupled BECs. The nonlinear energy is again e = 2.0. The coupling 
is K = 2.1, where both the m = ±1 and the m = ±2 are unstable. Until r = Tqsc ~ 11 we 
find very small oscillations of the relative particle difference which are of the order of 10~^. 
The decompositions into m modes shows that this is due to small oscillations taking place 
between the m = modes of the two rings (see also fig.[2]for r < 6). These are "ordinary" 
Josephson oscillations. Here no formation of currents, i.e ^L^^'^^ 7^ 0, takes place in either of 
the annuli. For r > 1 1 the situation changes dramatically. We find particle oscillations up to 
\N'^ — N^\/Ntot = 0.4 and a nonzero expectation values for l^L^^'^'^. The time Tosc of the onset 
of this regime depends on the magnitude of the initial seed A rough estimate for Tqsc 
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FIG. 3: (Color online) Angular momentum per particle in units of h and relative particle difference 
between the two annuli (e = 2.0 and k = 2.1). For r < 11 small oscillations in the particle 
difference in the order of 10"^ take place. At r = Tqsc ~ H (dashed vertical line) we observe 
the onset of oscillations with bigger amplitude which are accompanied by angular momentum 
Josephson oscillations. 

can be obtained from the equation A^o ~ = |«m(i = O)pexp(rrosc)- These oscillations 
in (Vi^'^'^ are due to the population of uj+ modes caused by the nonlinear evolution of the 
system. 

The oscillations spontaneously emerging for r > Tqsc are Josephson oscillations of the 
angular momentum. They break the chiral symmetry of the initial state's wave function 
where none of the rings carried a net angular momentum. 

Finally, we turn to the experimental realizability of the "ordinary" and angular momen- 
tum Josephson oscillations. For ^''Rb and a ring radius of po = 1-2 ^ifn the energy scale 
given by the length of the ring evaluates to -Eq ~ x 2nK and consequently we find a 



time scale of tq ~ 4 ms. For a radial oscillator length o 



oi e = 2.0 is achieved for a particle number A^o ~ 50 [16|]. The experimental feasibility of 



0.3 /xm a nonlinear energy 



building ring-shaped traps has been demonstrated recently 



[lO 



We therefore hope the 



results presented in this letter might stimulate further experiments. 

The required smallness of the ring traps forbids direct in-situ imaging of the BECs. 
Therefore one has to employ time-of-flight (TOF) imaging. We consider an experiment 
where initially two annular BECs are created in two uncoupled ring trap from one single 
BEC. Subsequently the barrier in z-direction is lowered such that a certain coupling strength 
K is established. The system then evolves at constant n for a certain time Tint after which 
the trap is switched off. The TOF image is then taken after free expansion of the cloud. 
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FIG. 4: (Color online) Simulated TOF images for several interaction times. The parameters 
correspond to those of the top panel of figure [2j The axis are in arbitrary units. For an interaction 
time of Tint = 0.0 only the m = mode contributes to the momentum distribution yielding an 
image of the squared zero-th Bessel function. At Tint = 10.7 the m = ±1 are occupied as well 
which gives rise to an angular modulation proportional to sin^((/)) of the TOF image. For later 
times several angular momentum modes are occupied. We show an example for Tint = 63.5. Here 
both annuli carry a net angular momentum. 



Here we assume that only the wave function of one annulus is imaged, i.e. the atoms in 
the second annulus have to be removed |17|. The TOF method yields an image of the 

I, uJ^ . This is equivalent to the squared modulus 



momentum distribution of the BEG 
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of the Fourier transform of the wave function \l/(r). Using the angular momentum mode 
decomposition am of the annulus which is to be probed and assuming that ap,az R we 
obtain 



vl/(k) oc Yl {-ir^'J\m\{kR)ame''^< - J] i-lY"'-'^/'j\m\{kR)ame 



imC, 



m=odd 



with Jn being the n-th Bessel function of the first kind, k = ^/c^ + ky and ( = aictan^ky/k^). 
Hence imaging |\l/(k) p allows the reconstruction of the am for a small number of contributing 
modes (for an example see fig.Hl). This allows an experimental study of the instability regions 
simply by a analyzing TOF images. 

In summary, in a system consisting of two ground state BECs in coupled rings the oc- 
cupation number of high angular momentum modes grows exponentially for well-defined 
coupling strengths. For small evolution times, a symmetric occupation of ±m modes takes 
place in each BEG accompanied by "ordinary" Josephson oscillations of the relative particle 
number. For later times, angular momentum Josephson oscillations spontaneously emerge. 
This novel type of Josephson oscillations breaks the initial chiral symmetry of the individual 
BEGs. 
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This can be achieved by first separating the rings and discarding one of them 
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